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A strongly interacting gas of two-electron fermions at an orbital Feshbach resonance 
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We report on the experimental observation of a strongly interacting gas of ultracold two-electron 
fermions with orbital degree of freedom and magnetically tunable interactions. This realization has 
been enabled by the demonstration of a novel kind of Feshbach resonance occurring in the scattering 
of two ^^^Yb atoms in different nuclear and electronic states. The strongly interacting regime at 
resonance is evidenced by the observation of anisotropic hydrodynamic expansion of the two-orbital 
Fermi gas. These results pave the way towards the realization of new quantum states of matter with 
strongly correlated fermions with orbital degree of freedom. 


Recent theoretical and experimental work has evidenced 
how quantum gases of two-electron atoms represent ex¬ 
traordinary systems for the development of a new gener¬ 
ation of quantum simulators of fermionic matter. They 
provide realizations of multicomponent fermionic gases 
with tunable number of nuclear spin states and SU(N)- 
symmetric interactions mi- In addition, they offer the 
unique feature of an additional electronic (orbital) degree 
of freedom, over which coherent quantum control can be 
achieved by means of the technology developed in the 
context of optical atomic clocks [4]. Recently, a strong 
spin-exchange interaction between ^^^Yb fermions in dif¬ 
ferent orbital states has been observed, paving the way 
to the realization of orbital quantum magnetism and of 
the Kondo lattice model urn- Despite these exciting 
perspectives, up to now two-electron atoms lacked the 
tunability of interactions that is provided by Feshbach 
resonances in the case of alkalis. Indeed, magnetic Fes¬ 
hbach resonances in ultracold gases of alkali atoms [ 7 ] 
have allowed breakthrough achievements, including un¬ 
precedented studies of strongly interacting fermions, with 
the demonstration of high-density molecular gases and 
the exploration of fermionic superfluidity at the BEC- 
BCS crossover [8]. A similar tunability for two-electron 
atoms would open totally new avenues, but the zero elec¬ 
tronic angular momentum in their ground state precludes 
the existence of accessible and useful magnetic Feshbach 
resonances. While optical Feshbach resonances in two- 
electron atoms have been proposed [9] and observed ex¬ 
perimentally [inilll], their actual implementation still 
suffers from severe intrinsic difficulties, such as heating 
and losses, preventing the observation of true many-body 
quantum physics, although novel promising schemes have 
been very recently investigated [12]. 

In this work we report on the first reaiization of a 
strongiy-interacting gas of two-eiectron fermionic atoms 
with orbitai degree of freedom and tunabie interactions. 
We take advantage of a recentiy-proposed orbital Fesh¬ 
bach resonance Ha affecting the scattering between ^^^Yb 
atoms in different eiectronic states. We observe a hy¬ 
drodynamic expansion of the Fermi gas in the strongiy 
interacting regime and use it to identify the resonance 



FIG. 1: Scheme of the ^^^Yb orbital Feshbach resonance 
affecting the scattering between one \g) — ^Sq atom and one 
|e) = ^Po atom in two different spin states (see text for details). 
The inset shows the behavior of the scattering length for 
Am = 5 as modeled in Ref. m, with the center of the 
resonance adjusted to match our experimental observations in 
Fig. 3 (the circle identifies the predicted position of the zero 
crossing). 


position. We also verify the predicted scaling of the reso¬ 
nance centers for different spin mixtures, arising from the 
SU(N) symmetry of ^^^Yb Fermi gases, and characterize 
the atom losses across the resonance, showing lifetimes 
that could allow for the investigation of the BEC-BCS 
crossover in two-orbital binary mixtures. 

The physical idea behind an orbital Feshbach resonance 
(OrbFR) is similar to that of a magnetic Feshbach res¬ 
onance, although with some distinctive differences. We 
consider two ^^^Yb atoms in two different electronic (or¬ 
bital) states 1^) = ^So (stable) and \e) = ^Pq (metastable, 
^ 20 s lifetime), and different nuclear spin states \f) 
and 1^). When the atoms are separated, interactions 
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between them are negligible and the relevant two-body 
eigenstates are \o) = f; e |) and |c) = that 

we name as open and closed collisional channels, respec¬ 
tively (see Fig. Q. The energy separation between the 
two channels is given by the differential Zeeman shift 
AgB = dgfiNAmB (where 6g is the differential Lande 
factor between \e) and |^), and /ijv is the nuclear magne¬ 
ton). In ^^^Yb A/i = 113 Hz/G x Am, where Am is the 
difference between the spin quantum numbers and m^. 
As the interatomic distance decreases, the appropriate 
basis for the description of the scattering is given by the 
orbital symmetric \eg^) = (|c) — |o))/\/2 and antisym¬ 
metric states \eg~) = (|c) + |o))/\/2, which are associated 
to two distinct molecular potentials, giving rise to two 
very different scattering lengths a+ and a“^, respectively 
0161. The difference between them determines an effec¬ 
tive coupling between open and closed channels, which 
becomes resonant when the differential Zeeman energy 
AjiB equals the binding energy Ec of the least bound state 
in the closed channel. Luckily, the orbital-symmetric po¬ 
tential features a shallow bound state with \Ec\ in the tens 
of kHz range, which corresponds to a large background 
scattering length a+. Thanks to this peculiar molecular 
potential of ^^^Yb, it is possible to use very convenient 
fields of the order of only tens of Gauss to magnetically 
tune the closed-channel bound state to be resonant with 
the open channel, despite the extremely weak sensitivity 
of two-electron atoms to external magnetic fields (about 
10^ times smaller than in alkalis). 

An ultracold gas of A ^ 6 x 10^ fermionic ^^^Yb atoms 
is prepared in a crossed optical dipole trap, operating at 
the magic wavelength = 759 nm in order not to shift 
the ultranarrow ^ ^Pq clock transition used for the 
preparation of the initial state. One of the trap beams 
is tightly focused, resulting in a cigar-shaped trap with 
angular frequencies uj = {ux^ojy^oUz) = 27r x (13,188,138) 
Hz. By proper spin manipulation protocols (relying on 
optical pumping and spin-selective removal of atoms by 
resonant excitation on the ^Pi transition), the 

gas is initially prepared in a balanced spin mixture of 
N/2 + N/2 ground-state atoms in nuclear spin states 
and m^, with different Am = — m^. The temperature 

of the Fermi gas is T 0.15 T;? (where Tp is the Fermi 
temperature) and the peak density is n 2.4 x 10^^ cm“^ 
per spin component. In order to excite the atoms to the |e) 
state, we first slowly turn on the intensity of a ID optical 
lattice at the magic wavelength with depth 30Ep (where 
pR = In?/2m}?^ is the recoil energy, with m the atomic 
mass). Along the direction of the lattice, we shine a 400 
/is resonant pulse of light on the 578 nm clock transition 
'5o ^ 3Po m, in such a way as to perform the optical 
excitation in the Lamb-Dicke regime [4]. This allows us 
to excite the atoms to the \e) state with an efficiency 
> 80% without imparting an optical momentum kick, 
which is important to maintain the system at equilibrium. 
The excitation is performed at a large magnetic field Be^c 
ranging from 72 to 167 G, in order to clearly resolve the 
Zeeman structure of the transition and excite only one 




FIG. 2: Anisotropic expansion of a strongly interacting ^^^Yb 
Fermi gas, initially prepared in the open channel of a spin 
mixture with Am = 5, at a magnetic field H = 41 G. a. 
The circles show the aspect ratio of the expanded \g) atomic 
cloud as a function of the time of flight. The long-dashed 
grey line shows the behavior expected for an ideal Fermi gas. 
The short-dashed black line indicates the unit limit above 
which the experimental aspect ratio gets inverted because of 
hydrodynamic expansion, b-d. False-color absorption images 
of the 1^) component of the atomic cloud after excitation in 
the open channel for increasing time of flight (b. 7 ms c. 16 
ms d. 28 ms). 


spin state. In this way we can selectively access the open 
or the closed channel (see Fig. [^. After the excitation, 
we adiabatically switch off the optical lattice in 100 ms to 
recover a 3D atomic gas. Then, we change the magnetic 
field intensity to the desired value B in about 2 ms. The 
trap is suddenly switched off and, after a time of flight 
tpoF^ the atoms remaining in the ground state \g) are 
imaged with a resonant pulse of light along i on the 
^ Pi transition. During the first 5 ms of the expansion 
the magnetic field is kept at the B value, thus allowing 
the atoms to release their interaction energy into kinetic 
energy. 

Fig.|2}i shows the evolution of the aspect ratio of the 
Fermi gas after the release from the trap as a function 
of the time of flight. The gas is prepared in the open 
channel of a spin mixture (m^ = — 5/2,m^ = +6/2) with 
Am = 5. The aspect ratio is defined as the ratio Ry/R^ of 
the expanded atomic cloud size along y to the size along 
X. In the case of a noninteracting Fermi gas (long-dashed 
lines) the expansion is ballistic, eventually resulting in 
a spherical shape and in an aspect ratio value of 1 for 
sufficiently long expansion times (much larger than the in- 
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verse trap frequencies). The experimental circles, showing 
the behavior of the interacting spin mixture at a mag¬ 
netic field 5 = 41 G, clearly show an inversion of the 
cloud shape from prolate to oblate, with an aspect ratio 
exceeding 1 after a Itof — 18 ms. Figs. -d show false- 
color absorption images of the atomic cloud for different 
times of flight, as specified in Fig. [^. The observation of 
the aspect ratio inversion is a hallmark of hydrodynamic 
expansion of the Fermi gas, which occurs in the regime of 
strong interactions, as observed for alkali fermionic gases 
close to magnetic Feshbach resonances m- In the hy¬ 
drodynamic limit the collisional rate F is larger than the 
geometric trapping frequency O = ^ causing 

a faster expansion along the tightly confined axis of the 
harmonic trap because of the larger density gradient. The 
observation of the aspect ratio inversion is an unambigu¬ 
ous signature of a collisional hydrodynamic regime and of 
strong interactions, which is a necessary condition for the 
onset of fermionic superfluidity pTlITS] . The quantitative 
value of the aspect ratio at resonance is lower than that 
predicted by the hydrodynamic equations of superfluids 
at resonant interactions (~ 2 at our maximum expansion 
time) [19]. This discrepancy has been already observed 
in m, and could be ascribed to the expanded gas not 
being fully in the hydrodynamic regime, possibly due to 
the narrow character of the OrbFR [21]. 

We observe that the asymptotic value of the aspect 
ratio in the open channel depends on the magnetic held 
value at which the expansion takes place, which is a strong 
evidence for the existence of the OrbFR predicted in Ref. 
m- Indeed, the anisotropy can be used as a tool to 
characterize the Feshbach resonance, as already studied 
with alkali atomic gases 120]. In Fig. we show the 
difference in the aspect ratio between the open channel 
and the closed channel for a spin mixture with Am = 5 as 
a function of the magnetic held at Itof = 28 ms. In the 
open channel a clear resonant behavior is observed, with a 
maximum that is located at a magnetic held 5 = (41 ± 1) 
G, signalling the enhancement of the elastic collisional 
rate at the Feshbach resonance. While we exclude possible 
confinement-induced shifts because of the low trapping 
frequencies in the 3D trap geometry, finite-energy effects 
could affect the resonance position since \Ec\ ^ h x 20 
kHz is only about one order of magnitude larger than the 
temperature and the Fermi energy kHz of 

the gas [13] (the latter corresponds to a magnetic held B = 
ksTp/Ajii 7 Gfor Am = 5, which might also contribute 
to the width of the observed features). Nevertheless, from 
the fitted value of the resonance center, using Eq. (13) of 
Ref. [13], we can extract afg -h a~g ^ 2100ao. However, 
we must note that this is only a rough estimate, given the 
narrow character of the observed resonance [7]. 

A more precise value would require further theoretical 
and experimental investigation, that will be postponed 
to future works. The value of a+ -h a~g also allows for 
an estimation of the magnetic held value at which the 
zero crossing of the scattering length occurs. Including 
it in Eq. (12) of Ref. [13], we obtain the behavior of the 



FIG. 3: Aspect ratio of the atomic cloud after a time of flight 
troF — 28 ms. a. Comparison between expansion in the open 
channel and expansion in the closed channel for a mixture 
with Am — 5. The solid line is a fit with a Lorentzian function. 
The aspect ratio background signal is less than unity because 
the expansion is not fully in the far-field limit, b Comparison 
between open-channel mixtures with different Am, plotted as 
a faction of a rescaled magnetic field B = BAm/5. Different 
colors refer to different spin states, as explained in the legend. 
The collapse of the different datasets onto the same curve is a 
verification of the OrbFR scaling law, which is in turn a direct 
consequence of the SU(N) symmetry of two-electron atoms. 
For these measurements we have used a different excitation 
scheme (the atoms are excited directly in the 3D trap without 
the applicaton of the optical lattice) and a different trap 
geometry d; = 27r x (22,181,139) Hz, which cause the peak 
aspect ratio to be different from that of the data in panel a. 
c. Resonance centers of the data shown in b (before rescaling 
the magnetic field) plotted versus Am. The line is a fit with 
the expected Am~^ behavior. 


scattering length for Am = 5, which is plotted in the inset 
of Fig. 12 The curve shows a zero crossing point around 
400 G, which lays outside the range of our experimental 
setup (170 G). 

We have repeated the measurements for smaller Am, 
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evidencing similar resonances at significantly larger mag¬ 
netic fields, as shown by the resonance centers plotted in 
Fig. The data are plotted in Fig. versus a rescaled 
magnetic field B = BAm/b. The different datasets clearly 
show a very similar dependence on 5, which is a distinc¬ 
tive feature of the OrbFR mechanism. As a matter of fact, 
the magnetic field at which the resonance is located is 
expected to scale as Am' m- This scaling law directly 
follows from the SU(N) invariance of scattering in ^^^Yb, 
which determines Ec to be independent of Am. As a con¬ 
sequence, the Zeeman energy dgii^AmB for which the 
coupling between the scattering channels is resonant, has 
to be the same for different Am, from which the scaling 
law B ^ Am~^ follows. The collapse of the experimental 
data onto the same curve is, therefore, a verification of 
the OrbFR scaling law, which is in turn a strong evidence 
for the SU(N) symmetry of two-electron atoms. 

For future studies of many-body physics in this novel 
experimental system, it is important to assess the effect 
of inelastic collisions at the Feshbach resonance. To this 
aim, we investigate the two-orbital Fermi gas, initially 
prepared in the open channel, for different holding times 
thoid in the trap. Fig. shows the measured aspect ratio 
as a function of the holding time at 5 = 60 G in the trap, 
while Fig. shows the number of atoms Ng remaining 
in the \g) state. From these measurements one can infer 
several important conclusions. First, we observe that 
the aspect ratio exhibits a slow decay (with 1/e lifetime 
r ^ 380 ms) towards the value (slightly smaller than 1) 
expected for a weakly interacting Fermi gas for the finite 
time-of-flight Itof = 28 ms. On the same timescale, Ng 
is decreasing towards a nonzero value. This behavior can 
be interpreted as the result of inelastic collisions between 
\e) and \g) atoms, which empty the ground state and 
cause the two-orbital Fermi gas to become progressively 
less interacting and to abandon the conditions for being 
collisionally hydrodynamic. The asymptotic nonzero value 
of Ng in Fig. can be interpreted as an excess of \g) 
atoms, caused by imperfections in the preparation of the 
initial state, which remain in the trap after the |e) atoms 
have been lost because of inelastic collisions. 

We note that the smooth decay of the aspect ratio in Fig. 
& is an evidence of the absence of shape excitations of the 
atomic cloud (which should have occurred on timescales of 
the order of the inverse trap frequencies), which confirms 
the adiabaticity of the excitation procedure. We have also 
verified that in the 3D geometry employed in the experi¬ 
ment the inter-orbital spin-exchange dynamics observed 
in Ref. [5] does not take place because of the small spin- 
exchange interaction energy Vex- the asymptotic closed 
and open channels are well defined by the differential 
Zeeman energy AgB ^ Vex and close to the Feshbach 
resonance we do not observe any significant repopulation 
of atoms in the \g state of the closed channel. 

Finally, we investigate the lifetime of the Fermi gas 
as a function of the magnetic field across the OrbFR, as 
illustrated in Fig. B- The lifetimes are estimated from 
a fit of the experimental data with a single exponential 
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FIG. 4: a. Evolution of the aspect ratio for a Am = 5 
cloud, initially prepared in the open channel, after different 
holding times in the trap at R = 60 G and a fixed time of flight 
troF = 28 ms. b. Number of atoms for the same experimental 
parameters as in a. The lines in a and b are exponential fits, 
c. Inverse 1/e lifetime of the Am — 5 open-channel mixture 
as a function of the magnetic field. 

function (see e.g. the line in Fig. |^), used to globally 
quantify the losses without differentiating among the pos¬ 
sible decay processes, which include atom+dimer and 
dimer+dimer inelastic collisions [2^ |23] . We observe a 
strongly asymmetric lifetime with respect to the center 
of the Feshbach resonance, with an increase in the loss 
rate on the BEG side, as expected by the activation of 
the aforementioned inelastic channels. The presence of 
multiple decay channels cause the single exponential fit 
not to be fully satisfactory in the description of the data, 
especially on the BEG side of the resonance, and further 
investigation is needed in order to ascertain the fundamen¬ 
tal processes at the basis of the atom losses. Nevertheless, 
the observed lifetime at the resonance seems to be rather 
long, with a 1/e lifetime of ~ 350 ms measured at 5 = 41 
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G, making this OrbFR a promising tool for future in¬ 
vestigations of many-body physics by using two-electron 
atoms. 

In conclusion, we have reported on the first realization 
of a gas of strongly-interacting two-orbital fermions with 
resonant interactions. This observation is the starting 
point for a whole new range of experimental investiga¬ 
tions, ranging from the investigation of the BEC-BCS 
crossover in an ultracold gas of fermions with orbital de¬ 


gree of freedom, including the realization of novel forms 
of topological superfluids with spin-orbit coupling m , to 
the investigation of two-orbital Hubbard models [25] with 
tunable interactions. 
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